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ABSTRACT 


The  principal  objective  of  this  report  is 
to  analyze  existing  information  concerning  the  char- 
acteristics of  atmospheric  turbulence  and  with  this  as 
a basis  to  arrive  at  some  representative  figures  on  the 
resulting  scattering  of  sound  to  be  expected.  This 
report  first  discusses  the  nature  of  such  turbulence, 
the  extent  and  causes  of  the  anisotropy  near  the  ground, 
the  variation  in  its  character  with  altitude,  and  the 
extent  to  which  approximations  of  isotropy  are  valid. 
Other  micrometeoro logical  parameters  are  considered  and 
their  effects  on  sound  propagation  in  comparison  to 

I 

I 

velocity  fluctuations  estimated.  Reference  is  made  to 
available  experimental  data  for  which  a reasonable 
spectral  function  for  the  turbulence  has  been  construc- 
ted. This  function  is  applied  to  a scattering  cross 
section  expression  (8)  to  obtain  numerical  results 
for  the  scattering.  These  are  compared  with  the  few 
available  experimental  data. 


1.  Introduction 


Of  the  manifold  meteorological  disturbances 
that  affect  the  transrai ssion  of  sound  through  the  lower-  1 
most  layers  of  the  atmosphere,  the  turbulent  motion  of  the 
air  has  been  recognized  for  some  time  as  playing  a role  of 
utmost  significance.  The  characteristically  severe  fluc- 
tuations in  received  intensity  of  sound  which  results  from 
it  have  been  widely  commented  upon  by  numerous  observers 
and  some  qualitative  results  obtained.  Thus,  for  example, 
it  is  agreed  that  these  fluctuations  increase  with  the  in- 
tensity of  the  turbulence,  the  frequency  of  the  sound,  and 
the  separation  between  source  and  the  receiver  (1,  2).  Yet 
quantitative  records  of  such  phenomena  are  singularly 
limited  in  number  and  detail.  In  the  recent  literature, 
papers  by  Sieg  (3),  Krasilnikov  (4)>  and  Richardson  (jj>), 
and  by  Schilling  and  his  associates  on  ultrasonics  propaga- 
tion (6),  stand  out  as  isolated  examples. 

The  dearth  of  concrete  data  can  hardly  be  ascribed 
to  the  lack  of  important  applications  since,  apart  from  its  1 
intrinsic  interest,  such  information  is  of  value  in  studies 
i of  military  sound  ranging  and  aircraft  noise  propagation. 

The  reason  is  to  be  found,  as  a reading  of  the  literature 
readily  makes  apparent,  in  the  difficult  experimental  pro- 
cedures encountered  in  micrometeorological  research.  Indeed, 
as  Schilling,  et  al  (7)  point  out,  any  inquiry  into  the 
correlation  between  atmospheric  sound  propagation  and  local 
micrometeorological  conditions  immediately  draws  the  lnvesti- 
i gator  into  the  diverse  and  involved  problems  of  the  micro- 
meteorologist.  The  instrumental  demands  are  of  the  most 
exacting  nature,  e.g.,  the  sensory  elements  employed  (thermo- 
couples, anemometers,  etc.)  must  be  capable  of  recording 
accurately  fluctuations  whose  periods  extend  from  milliseconds 
to  many  minutes.  Careful  interpretation  of  results  is  called 
for  since  instruments  which  measure  the  same  parameter  may, 
through  differences  in  characteristics,  yield  data  which  are 
respectively  difficult  to  compare. 

The  situation  described  above  makes  it  difficult 
to  prepare  a set  of  consistent  data  in  the  form  of  tables  or 
graphs  which  would  allow  quantitative  predictions  to  be  made 
on  the  amount  of  turbulent  scattering.  The  solution  obviously 
lies  in  improved  instrumentation  and  further  experimental  work 


Theoretical  expressions  for  the  expected  scattering  in 
a turbulent  medium  have  been  derived  by  Kraichnan  (8,  9)* 
However,  the  application  of  these  equations  to  atmospheric 
turbulence  likewise  requires  definite  information  concern- 
ing its  spectrum,  and  the  above  mentioned  difficulty  is 
encountered  again. 

The  principal  objective  of  this  report  is  to 
analyze  as  much  information  concerning  the  characteristics 
of  atmospheric  turbulence  as  is  presently  available  and 
with  this  as  a basis  to  arrive  at  some  representative 
figures  on  the  scattering  to  be  expected.  Consequently 
this  report  will  be  divided  into  two  parts,  of  which  the 
first  is  a review  of  the  available  micrometeorological 
data  pertinent  to  the  study  of  sound  propagation  in  open 
air.  Such  data  will  be  derived  solely  from  the  available 
existing  literature,  since  no  new  experimental  work  is  to 
be  reported  upon  in  this  paper.  The  second  part  consists 
of  calculations  in  which  expressions  derived  in  the  the-  | 
oretical  papers  by  Kraichnan  (8,  9)  will  be  used  in  con- 
junction with  the  above-mentioned  data  to  indicate  the 
scattering  to  be  expected  for  the  case  of  a medium  iso- 
tropically turbulent,  which  represents  conditions  to  be 
found  in  the  upper  regions  of  the  atmospheric  boundary 
layer  and,  for  a certain  range  of  eddy  sizes,  conditions 
found  close  to  the  ground.  In  view  of  the  inexact  charac- 
ter of  the  data  employed,  the  figures  obtained  will  apply 
strictly  only  to  the  particular  case  ftonsidered,  although 
in  order  of  magnitude,  they  will  be  correct  for  a wider 
range  of  data. 

Before  proceeding  further,  mention  should  be 
made  of  factors  other  than  turbulence  which  contribute 
to  the  attenuation  of  airborne  sound.  The  absorptive 
effect  of  humidity  and  the  refractive  bending  of  wave 
fronts  due  to  the  vertical  stratification  of  wind  vel- 
ocity and  temperature  are  well  known  and  have  been 
reported  upon  in  a number  of  papers  (2,  10,  11,  12).  The 
latter  two  are  responsible  for  the  long  period  fluctua- 
tions which,  depending  upon  the  sign  of  the  gradients 
involved,  either  cause  or  prevent  the  formation  of 
"shadow  zones".  It  is  evident  that  the  short  period 
fluctuations  on  the  order  of  seconds,  which  are  the  concern 
of  this  paper,  do  not  involve  these  large  scale  vertical 
gradients,  but  only  the  high  frequency  departures  there- 
from. In  addition  to  the  attenuation  due  to  water  vapor 
In  the  air,  the  effects  of  viscosity  and  spherical  diverg- 
ence are  also  present.  Though  the  first  two  show  a 
rising  characteristic  with  frequency,  all  three  are  of 
negligible  importance  here.  Ground  attenuation  as  a 
factor  is  also  excluded  from  our  consideration  (12,  13). 
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2.  Notation 


, In  what  follows,  the  customary  notation  uc(i  =\,x, 3) 

as  denoting  instantaneous  wind  velocities  along  the  x,  y, 

, and  z directions,  respectively,  will  be  used.  Similarly, 

, will  refer  to  the  three  space  coordinates.  On 

i the  assumption  that  most  eddies  are  horizontal,  which  is 
frequently  the  case,  if  the  direction  of  the  mean  wind 
, velocity  is  chosen  in  the  x-direction,  cross-wind  will  be 
, along  the  X-axis  and  vertically  upward  will  be  the  *$axis. 

The  mean  wind  velocity  is  defined  by  the  expression  ir^T'* (ru, 

, where  T is  a suitably  long  time  interval  over  which  the  ° ’ 
i average  is  taken.  (It  is  acknowledged  that  the  specifica-  1 
, tion  of  such  a mean  value  involves  an  arbitrary  choice  of 
|T.  To  avoid  a lengthy  discussion,  we  shall  adhere  to  cus- 
i tomary  practice,  namely,  the  choice  of  T depends  upon  the 
periods  of  the  fluctuations  of  interest.)  According  to  the 
above  orientation  of  the  coordinate  axes,  ul  and  are  zero. 

The  instantaneous  eddy  velocities  u'.  are  defined  as  follows: 

< is  the  difference  between  the  instantaneous  wind  velocity 
u,  and  the  mean  wind  velocity  u,:  uf ^ ^ . 

The  parameter  which  is  correlated  to  the  scattering 
of  sound  waves  is  the  "gustiness"  or  "intensity  of  turbulence. 
This  is  identified  by  the  symbol  G,  and  is  defined  by  ■&&)'*'/ 
This  is  the  root  mean  square  eddy  velocity  divided  by  the  mean 
velocity. 

3.  Turbulence  in  the  Atmosphere 
3.1  Introduction 


Atmospheric  turbulence  is  essentially  a large 
scale  boundary  layer  problem  in  which,  in  contrast  to 
the  boundary  layer  depths  encountered  on  the  laboratory 
scale,  on  the  order  of  centimeters,  the  depth  of  the 
,layer  is  many  hundreds  of  meters.  This  gives  it  the 
advantage  of  instrumental  accessibility;  on  the  other 
hand,  the  bounding  surface  displays  irregularities  never 
found  in,  say,  a wind  tunnel,  a circumstance  which 
greatly  complicates  the  nature  of  the  resulting  turbul- 
ence. The  main  flow,  essentially  laminar,  consists  of 
the  freely  flowing  winds  In  the  atmosphere  above  the 
boundary  layer,  which  are  driven  by  pressure  gradients 
and  the  Coriolis  force.  Upon  coming  into  contact  with 
the  surface  of  the  earth,  these  generate,  through 
friction,  turbulent  eddies,  and  set  up  a vertical  strati- 
fication of  air  motion  in  the  lower  atmosphere.  These 
eddies  are  not  to  be  thought  of  as  simply  circular 
vortices  or  waves  but  rather  as  fluctuations  in  the  mean 
wind  velocity,  caused  by  moving  masses  of  air  with  non- 
zero vorticity.  As  a series  of  similar  eddies  move  past 
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the  observer,  a periodicity  in  wind  velocity  and 
wind  direction  frequently  will  be  observed.  An  es- 
timate of  the  "eddy  size'*  can  be  made  if  the  period 
and  the  mean  wind  velocity  are  known  (Figure  1). 

* Usually,  however,  the  velocity  trace  of  a turbulent 
flow  is  quite  irregular  (Figure  2)  and  only  rough 
estimates  of  size  and  period  can  be  made,  and  for 

* purposes  of  analysis,  a Fourier  decomposition  is 
generally  resorted  to  (Section  3) * 

The  atmospheric  boundary  layer  can  be 
divided  into  two  regions.  The  layer  lowermost,  the 
surface  boundary  layer,  extends  to  a height  of  about 
100  meters.  Being  immediately  adjacent  to  the  ground 
(or  ocean),  large  inequalities  are  set  up  in  the  vel- 
ocity fluctuation  components  along  the  three  coordinate 
axes.  Such  anistropy  of  turbulent  motion  is  quite 
characteristic  of  this  layer.  However,  among  those 
eddies  whose  size  is  small  enough,  of  order  x3/u(, 
isotropy  may  locally  exist.  That  is,  the  mean  square 
velocity  fluctuations  are  independent  of  the  orientation 
of  the  coordinate  axes:  h',’'  =«*.£  - . This  is  consistent 

with  the  observation  that  the  eddy  size  increases  with 
height  so  that  isotropy  of  eddy  motion  extends  to  pro- 
gressively larger  eddy  sizes  as  the  elevation  is  increased. 
Above  this  layer,  extending  to  a height  of  about  1 km,  is 

* the  planetary  boundary  layer  which  shows  little  anistropy, 
the  velocity  fluctuations  being  fairly  independent  of 
direction. 


It  is  evident  from  the  above  that  the  study  of 
acoustic  propagation  involves  the  tracing  of  ray  paths 
through  either  the  surface  layer  alone,  in  which  case 
the  medium  is  predominantly  anisotropic,  or  through  both 
layers,  the  choice  depending  on  the  particular  refractive 
conditions  prevailing.  In  the  interests  of  simplicity, 
however,  we  shall  consider  only  the  two  pure  cases  of 
complete  isotropy  and  anisotropy  in  this  report. 

The  specification  of  conditions  for  isotropy 
is  quite  simple;  consequently  the  main  body  of  micro- 
meteorological  data  to  be  considered  will  pertain ,to  the 
departure  from  Isotropy  which  is  to  be  found  in  the 
surface  boundary  layer. 

The  random  velocity  fluctuations  u- (?,-£)  which 
characterize  turbulent  flow  in  the  atmosphere  are  cus- 
tomarily assumed  to  arise  from  contributions  to  the  flow 
by  eddies  of  different  sizes.  Accordingly  a Fourier 
decomposition  of  is  resorted  to: 
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where  is  the  wave  number  vector  belonging  to 

the  length  £ , and  v-*  is  the  Fourier  amplitude*  can 
be  considered  to  be  a kind  of  eddy  size. 

^ The^  values  of  u'  at  two  separate  points  in 
space:  x and  *'  , will  define  a correlation  temsor : 

<r) 

In  the  event  the  turbulence^  is  isotropic,  tijwill  be  a 
function  of  the  separation  A - > ~'i'  only: 

tCx)  ~ ( k'tl ) t)  !*.'( £+  x 

The  Fourier  transform  of  4\>is  given  by: 

(x*)“>l  \ ^ 43* 

- 

and,  under  isotropic  conditions,  the  energy  density  of 
turbulent  fluctuations  that  lie  within  the  limits  k and 
k+<ifc  is: 

rr 

£(k)  - (*-*)  J J (k>  k*-  s,v,  e-  tie  dcf, 

o O 

An  idea  of  the  average  eddy  size  is  obtainable  by  in- 
tegrating ±(\)  as  a function  of  \ . From  the  definition 
of  the  correlation  function,  I when  V-o  and  monotonically 
decreases  to  zero  for  some  large  value  of  \ . The  more 
rapidly  +00  approaches  zero,  the  smaller  the  eddy  size. 
Therefore  we  shall  define  the  average  eddy  size  as 

l - J 4(»j  as 


3.2  Isotropic  Turbulence 

In  the  study  of  Isotropic  turbulence.  It  has 
been  found  convenient  to  look  for  situations  where 
measurements  are  possible  under  controlled  conditions. 
Advantage  is  taken  of  the  fact  that  the  assumption  of 
isotropy  reduces  the  significant  parameters  to  two: 
gustiness  and  scale.  One  such  opportunity  exists  in- 
side wind  tunnels  where  turbulent  flow  is  produced  by 
means  of  a fine  mesh  grid.  At  distances  * on  the  order 
Of  */*>  ~ iol  where  M is  the  mesh  size,  the  turbulence  decays 
to  an  approximately  isotropic  state  (lij.).  It  must  be 
kept  in  mind,  however,  that  conditions  found  in  a wind 
tunnel  differ  from  ground  conditions  in  that  there  is 
no  thermal  stratification  and  that  an  upper  limit  to  the 
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eddy  size,  which  is  determined  by  the  mesh  size  U3ed, 
exists.  This  tends  to  make  the  gustiness  observed  in 
wind  tunnels  smaller  by  a factor  of  five  or  more  than 
that  actually  found  close  to  the  ground  (15).  In  spite 
: of  this,  the  ease  with  which  the  above-mentioned 
parameters  can  be  varied  makes  the  wind  tunnel  an  in- 
dispensable source  of  information  in  turbulence  study. 

The  isotropic  state  as  it  exists  in  the  wind 
tunnel  is  closely  represented  by  the  3-dimensional 
Heisenberg  function 

= tV  (ulH1)"’ 


where  L is  the  scale  of  turbulence  defined  previously. 

It  is  easily  seen  that  for  small  wave  numbers 
and  for  k large,  k»L~'  ;e*kz(Fig.  3). 

3*3  Effect  of  the  Ground 

The  turbulent  motion  of  the  air  which  directly 
alters  the  character  of  sound  propagation  is  found  to  be 
most  intense  within  the  first  few  meters  of  the  ground. 

This  is  a consequence  of  the  dominant  role  played  by  the 
ground  in  the  formation  of  eddies  (a)  the  ground  is  ex- 
posed each  24  hours  to  a heating- cooling  cycle,  the 
periodic  temperature  variations  of  which  produce  alterna- 
ting positive  and  negative  temperature  gradients  immediately 
adjacent  to  the  surface;  (b)  irregularities  in  terrain 
features,  by  applying  frictional  drag  to  the  prevailing 
air  flow,  generate  turbulent  eddies  whose  mean  size  is 
characteristic  of  the  terrain.  A shear  flow  is  thus 
created  with  strong  components  parallel  to  the  surface 
contours. 


(a)  Thermal  Effects:  The  negative  vertical 

temperature  gradient  Vr/n *.o{  normal  lapse  rate)  which 
occurs  during  the  day  results  In  a convective  upward 
transport  of  air  which  is  superposed  upon  the  prevailing 
wind  and  is  maximum  near  the  ground  about  noon.  Accom- 
panying this  upward  flow  Is  usually  an  inflow  of  denser 
cool  air  which,  upon  striking  the  ground,  produces  short- 
period  frictional  eddies.  Subsequently,  the  cool  air  is 
warmed  by  the  ground  and  itself  rises.  It  is  apparent 
that  this  instability  in  density  produces  a large  scale 
circulatory  motion  of  the  air  in  a vertical  plane.  This 
motion  takes  on  a cell-like  aspect,  which  has  been  studied 
by  Bfcnard,  and  by  Rayleigh  and  others,  who  conclude  that 
the  characteristic  size  of  the  eddies  thus  formed  is  of 
the  order  of  1 km  in  the  direction  of  the  wind  and  about 
half  as  much  across.  It  is  to  be  noted  that  the  latter 
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figure  exceeds  by  many  times  the  wavelength  of  sounds 
normally  encountered  and  hence  that  this  type  of  eddy 
is  not  likely  to  be  important  in  producing  scattering. 
However,  long  period  fluctuations  in  intensity  commonly 
known  as  fading  are  ascribed  to  such  large  eddies.  It 
may  thus  be  concluded  that  the  generation  of  frictional 
eddies  by  Benard  cells  is,  for  this  study,  the  most 
important  result  of  the  diurnal  heating  cycle.  Paren- 
thetically it  should  be  added  that  the  thermal  stratifi- 
cation which  the  lapse  rate  represents  is  itself 
responsible  for  refractive  effects  in  which  distortion 
of  the  acoustic  wave  front  occurs,  especially  in  the 
vertical  direction.  As  the  lapse  rate  decreases  with 
approach  of  night,  there  is  a corresponding  decrease  , 

in  the  gustiness.  The  transition  from  a normal  lapse 
rate  to  inversion  *° that  occurs  as  a result  of  rapid 
cooling  of  the  ground  causes  a reversal  of  the  density  , 
gradient  and  stabilizes  the  air  motion.  The  large 
eddies  referred  to  die  out  completely,  with  an  accom- 
panying reduction  in  the  intensity  of  the  frictional 
eddies.  That  the  latter  do  not  vanish  altogether  is  due 
to  the  persistence  of  the  mean  flow.  However,  under 
extreme  inversional  conditions,  it  is  found  that  even 
the  frictional  eddies  are  damped  out.  The  fact  that  at 
night  sounds  are  heard  more  clearly  is  due  (aside  from 
the  general  reduction  in  ambient  noises)  partly  to  this 
decrease  in  gustiness,  although  the  more  favorable  re- 
fractive properties  of  the  night  air,  which  inhibits 
the  formation  of  "shadow  zones"  is  the  predominant  factor. 

(b)  Terrain  Effects:  The  nature  of  the  earth 

surface  affects  sound  propagation  by  its  role  in  the  pro- 
duction of  frictional  eddies  mentioned  above,  and  by  the 
absorptive  and  reflective  role  of  the  ground  covering. 

The  fluctuations  In  intensity  and  phase  of  interest, 
however,  are  due  to  the  first  factor.  Since  these 
frictional  eddies  arise  from  the  flow  of  air  over  ob- 
stacles such  as  grass,  trees,  and  buildings,  it  is  natural 
to  expect  that  vorticity  not  be  oriented  at  random,  but 
for  the  fluctuations  In  the  direction  of  the  wind  to  be 
greater  than  in  other  directions.  This  is  shown  in  the 
next  section  to  be  the  case.  It  is  not  easy  to  correlate 
the  properties  of  the  ground  with  the  quantities  commonly 
associated  with  turbulence,  but  according  to  Sutton  ( 16 ) , 
the  gustiness  should  be  related  to  the  drag  coefficient 
Cd.  The  accompanying  table  shows  some  representative 
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values  for  C^.  It  can  be  seen  that  Cd  is  proportional 
to  what  one  might  consider  the  roughness  of  the  surface. 


Nature  of  Surface 

Cd 

Very  smooth  (mud  flats,  ice) 

.002 

Grass  to  1 cm 

.005 

Thick  grass  to  10  cm  high 

.016 

Thick  grass  to  50  cm  high 

.032 

Sea 

.002 

Neutral  equilibrium: 

it,  - Too  (w/ju.. 

The  role  that  the  ocean  plays  in  contributing 
to  turbulent  conditions  above  its  surface  Is  quite  com- 
plicated since  the  roughness  of  the  ocean  is  Itself  a 
function  of  the  wind  velocity. 

The  drag  coefficient  Cd  for  the  ocean  given 
i in  the  table  above  (Ct-.ooi)  applies  only  under  the  con- 
1 ditions  stated  of  neutral  equilibrium  and  wind  speed  of 
1 5m  sec'.1  Observations  seem  to  indicate  that  Cd  is  pro- 
1 portional  to  the  wind  speed  and  that  for  conditions  other 
1 than  neutral  equilibrium,  other  values  of  Cd  must  be 
1 employed.  The  percentage  turbulence  over  the  ocean  seems 
to  either  rise  or  remain  constant  with  height,  which 
1 indicates  the  persistence  of  turbulence  to  heights  not 
found  on  land.  There  is  evidence  (25)  to  Indicate  that 
the  height  to  which  the  turbulence  extends  above  the 
ocean  includes  the  tropopause,  and  perhaps  into  the 
stratosphere.  In  addition  to  the  roughness  due  to  shear 
motion  of  the  air,  the  part  played  by  the  turbulent 
motion  of  the  ocean  itself  in  causing  air  turbulence  is 
not  known.  There  is  evidently  in  this  case  very  little 
information  with  which  to  work. 

3.k  Anisotropy  of  Eddy  Velocities 

It  is  a well  established  fact  that  close  to 
the  ground,  fluctuations  in  the  vertical  direction  are 
not  as  great  as  in  the  horizontal.  This  difference, 
however,  becomes  less  with  increasing  height  (18).  Such 
behavior  is  usually  expressed  meteorologically  by  the 
specification  of  the  Reynolds  shear  stresses:  «-/•»*< 

where/0  is  the  density.  These  stresses  are  indicative 
of  momentum  transport  by  the  velocity  fluctuations 
across  the  surfaces  defined  by  i and  j.  A non-zero 
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value  of  implies  a correlation  between  u(  and  u'} 
at  that  point,  and  therefore  a preferential  transfer 
of  momentum.  In  particular,  t(5=  -p u5  called  the 
horizontal  shear  stress,  measures  the  momentum  trans- 
port across  any  plane  parallel  to  the  ground.  Accord- 
ingly, the  existence  of  a non-zero  value  ofTr^in  a 
region,  indicating  correlation  between  horizontal  and 
vertical  eddy  motion,  is  a condition  for  the  presence 
of  anisotropy.  This  situation  has  been  observed  in, 
and  actually  defines  the  extent  of,  the  atmospheric 
boundary  layer. 

Observations  on  the  relative  difference  , 

between  the  various  eddy  components  show  that  invari- 
ably u'j  exceeds  u',  . Scrase  (19)  reported  that  under 
conditions  of  small  temperature  gradient  < was  larger 
than  by  $0%  at  a height  of  2 meters  with  a gradual 
decrease  to  20  meters,  after  which  no  asymmetry  was 
found.  This  decrease  has  been  questioned  by  several 
investigators  (20)  but  no  satisfactory  explanation  has 
been  forthcoming.  Best  (21)  found  on  the  other  hand, 
at  a height  of  200  meters,  that  G2  exceeded  G3  by  about 
80$.  This  indicates  that  large  variations  in  aniso- 
tropic conditions  may  be  expected  under  varying  terrain 
conditions.  The  importance  attached  to  the  ratio  ui:u$ 
lies  in  its  effect  upon  the  vertical  and  horizontal  phase 
relationships  in  received  sound  which  are  especially 
significant  in  acoustic  direction  ranging.  It  is  re- 
ported by  Richardson  (5)  that  phase  fluctuations  were 
noticeably  greater  when  a source  and  receiver  were  set 
6 ft.  apart  vertically  than  when  set  12  ft.  apart 
horizontally.  A similar  effect  is  also  reported  by 
Sieg  (3).  „ 


A representative  set  of  data  illustrating  the 
anisotropy  is  given  in  the  accompanying  table  (22). 


MEAN;  EDDY 

VELOCITIES  AND 

SHEARING  STRESS 

^Period  (EST) 

V 

*3 

u^ru^u^ 

1440-1447 

140 

88 

32 

1.0:0.63:0.37 

4.5 

1522-1529 

122 

83 

50 

1.0:0.68:0.41 

4.8 

1705-1712 

103 

70 

51 

1.0:0.68:0.50 

4.3 

1729-1736 

92 

60 

56 

1.0:0.65:0.61 

5.0 

2110-2117 

45 

25 

15 

1.0:0.56:0.33 

0.2 

2150-2157 

63 

36 

23 

1.0:0.57:0.37 

1.1 

Taken  at  height 

of  2.3 

meters. 

It  can  be  seen  that  in  general  m.',  > and 

that  T(1  decreases  with  the  change  from  normal  lapse 
rate  to  inversion  conditions.  Data  taken  by  other  in- 
vestigators show  a similar  range  of  values  for  'b,*,  the 
differences  being  ascribed  to  mean  wind  speed,  rough- 
ness of  terrain,  the  height  of  observation,  thermal 
stability  of  the  atmosphere,  and  not  least,  the  in- 
strumentation employed. 

As  might  be  expected,  the  probability 
distribution  of  the  eddy  velocity  amplitudes  is  Gaussian 
(21,  23) • Each  component  roughly  follows  an  equation 
of  the  form: 

The  gustiness,  defined  previously,  is  approx- 
imately independent  ofG,.  It  is  shown  (22,  23)  that  the  , 
velocity  fluctuations  increase  in  proportion  to  the  mean  , 
wind,  a dependence  which  seems  to  persist  over  all 
ranges  of  the  lapse  rate.  Best  (21)  has  found  that  G^ 
is  approximately  constant  for  2.5**j^  200  cm,  and  is 
fairly  independent  of  temperature  in  the  range  2.5  * < 

10  cm,  with  however  a decrease  as  conditions  change  from 
lapse  to  inversion.  Theory  (2i|)  indicates,  on  the  other 
hand,  that  G,  increases  close  to  the  ground,  as  conditions 
tend  toward  inversion.  This  seems  to  be  supported  by 
Cramer's  data  (22)»  There  i3  evidently  more  data  necessary 
before  a definite  conclusion  is  reached.  Representative 
values  of  Gi  given  in  the  literature  range  from  0.1  to  O.I4.. 

3.5  Eddy  Sizes 

Observations  by  Scrase  (19)  and  Panofsky  (25) 
indicate  that  the  average  eddy  size  generally  increases 
with  altitude,  which  is  reasonable  in  view  of  the  decrease 
in  frictional  agents.  The  scale  of  turbulence  L which  is 
a rough  measure  of  the  average  eddy  size,  owes  its  im- 
portance to  the  fact  that  the  wave  front  of  a compressional 
wave  is  distorted  when  it  encounters  a shear  flow,  the 
magnitude  of  the  resulting  scattering  being  determined  by 
the  size  of  the  flow.  We  shall  consider  the  longitudinal 
scale  L, measured  in  the  mean  direction  of  wind  and  the 
transverse  scale  Lxmeasured  in  a crosswind  direction.  It 
has  been  observed  that  L,is  often  longer  than  Lx(22).  As 
a matter  of  fact,  there  seem  to  exist  daily  fluctuations 
In  both  scales  in  which  L,is  larger  than  L-v.  during  the 
early  afternoon  and  at  night,  while  during  the  late 
afternoon  they  are  approximately  equal.  At  night  both 
show  a decrease  in  size  (22).  The  representative  values 
of  L, taken  at  a height  of  2.5  meters  range  from  about  6m 
at  night  to  about  8m  in  late  afternoon  whereas  these  of 
Lx  ranges  from  about  7m  during  the  same  period.  This  is 
indicative  of  an  average  eddy  size  of  less  than  10  meters 
and  is  in  agreement  with  the  scale  found  by  Schilling 
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et  al  (6).  with  a mean  daytime  wind  velocity  of  £ m/sec. , 
TiEis  would  correspond  to  a mean  eddy  period  of  1 to  2 sec. 
The  reduction  in  mean  wind  velocity  at  night  to  about 
1 m/sec  increases  the  eddy  period  to  about  6 sec.  These 
values  agree  fairly  well  with  those  found  by  Scrase  (19) 
and  others. 

It  is  to  be  noted  that  measurements  of  scale  as 
calculated  above  depend  very  greatly  upon  the  time  inter- 
val over  which  the  observations  are  taken.  This  of  course 
is  due  to  the  increasing  range  of  eddy  sizes  included.  It 
has  been  found  that  as  the  time  increases,  the  correlation 
coefficients  rise  asymptotically  to  a constant  value.  A 
representative  value  for  the  time  in  which  such  asymptotic 
values  are  reached  by  the  longitudinal  and  transverse 
coefficients  is  about  100  seconds  (22). 


Calculations 

In  the  paper  by  Kraichnan  previously  referred 
to  (8),  an  explicit  expression  for  the  averaged  total 
cross  section  for  scattering  for  isotropic  turbulence 

is  given.  According  to  this  formula: 

<<r>  * (5. 16)  (8) 

I c o 

where  <r  is  the  total  scattering  cross  section  per  cm-', 

1 is  r.-rr  times  the  incident  frequency,  c is  the  velocity 
of  sound,  x is  sin  e is  the  angle  between  the  incid- 

ent wave  and  the  scattered  wave,  and  E(2«o«x/c)  is  the 
spectrum  function  for  the  kinetic  energy  per  unit  mass. 

The  three  dimensional  Heisenberg  spectral  density  for 
isotropic  turbulence  previously  given  is: 

EOO  - Uk */(<+ 

where  L is  the  scale  of  the  turbulence  and  k=  77  is  the 
wave  number  of  the  flow.  In  order  to  employ  this  function 
in  the  expression  for  <<r>  , it  is  necessary  to  form  the 

product  where  is  the  mean  square  of  the 

velocity  fluctuation  components. 

Since  K-^^the  quantity  is  given  by: 

<EClrjr')>=  LV  yyo-  lV)’  . 
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The  complete  expression  for  <cr>is  then: 


<«->  = 


I 4*  — *r—  * 

c 


•r 


where  G(x)  is  the  quantity  plotted  in  Figure  2, 
(8),  and  is  defined  as: 


reference 


QU)  = IxSO-i+fO-x'-)''- 

Since  the  principal  scattering  occurs  for  eddy  sizes  on  the 
order  of  the  incident  wavelength  , we  may  set  X#~L  so 

that  s X ~ . In  which  event,  <<r>  becomes 

*«r>  = d* 

% (l-*-  »C>-rr" **■)'’ 


It  is  apparent  that  this  expression  indicates  a 
'direct  dependence  of  the  scattering  upon  the  incident 
'frequency  fQ  and  the  gustiness,  in  agreement  with  exper- 
imental evidence.  It  will  be  recalled  that  the  gustiness 
was  previously  defined  as  the  root  mean  square  velocity 
fluctuation  (u.il)Vl  divided  by  the  mean  velocity. 


Numerical  .integration  of  the  integral  yields  the 
value  of  3.5  * 10“5.  With  typical  values  of  fp  - 100 
cycles,  u7‘  =■  1.5  * lo4  cm^sec-^  (23) » and  c = 3*4-  *10^  cm/sec, 
we  get  for<<r>}1.6  * 10“°  cm“l.  Emplbying  the  expression 
I = I exp  (-ax)  whore  I is  the  intensity  at  distance  x, 
and  X®  is  the  initial  intensity,  calculation  shows  the  loss 
per  foot  to  be  0.21  * 10”3  db  at  100  cycles. 


The  table  below  presents  results  from  sin^ilar 
calculations  for  several  frequencies.  For  comparison, 
experimental  results  recently  obtained  by  Delsasso  (27)  in 
the  high  Sierras  under  varying  atmospheric  conditions  are 
shown  alongside.  The  elevations  at  which  these  data  were 
taken,  about  500  m.  maximum,  are  high  enough  to  assume 
that  isotropy  exists  for  the  eddy  sizes  under  consideration. 
These  values  are  not  strictly  comparable,  however,  since 
the  attenuation  due  to  diffraction  and  humidity  is  also 
included. 


frequency 

calculated 

Results 

observed  (Delsasso, 2?) 

100  ^ 

0.21  * 10*3  .012  to  1.12  * 10"3 

25  o 

0.53 

.36  1.47 

500 

1.05 

1.06  2.64 

1000 

2.10 

.83  2.56 

It  is  apparent  that  the  expression  from  which  the  above  re- 
sults were  calculated  yields  results  of  the  correct  order 
of  magnitude.  ' - 
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